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Abstract 

Properties of the phase space of the standard map with memory are investigated. This map was obtained from a 
kicked fractional differential equation. Depending on the value of the map parameter and the fractional order of the 
derivative in the original differential equation, this nonlinear dynamical system demonstrates attractors (fixed points, 
stable periodic trajectories, slow converging and slow diverging trajectories, ballistic trajectories, and fractal-like 
structures) and/or chaotic trajectories. At least one type of fractal-like sticky attractors was observed. 
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1. Introduction 

The standard map (SM) can be derived from the dif- 
ferential equation describing kicked rotator. The de- 
scription of many physical systems and effects (Fermi 
acceleration, comet dynamics, etc.) can be reduced to 
the studying of the SM |[l| . The SM provides the sim- 
plest model of the universal generic area preserving map 
and it is one of the most widely studied maps. The 
topics examined include fixed points, elementary struc- 
tures of islands and a chaotic sea, and fractional kinetics 

It was recently realized that many physical systems, 
including systems of oscillators with long range inter- 
action [4, JJ, non-Markovian systems with memory (|01 
Ch.lO, HMHEIilll), fractal media etc., can be 
described by the fractional differential equations (FDE) 
[H, As with the usual differential equations, 
the reduction of FDEs to the corresponding maps can 
provide a valuable tool for the analysis of the proper- 
ties of the original systems. As in the case of the SM, 
the fractional standard map (FSM), derived in lll5ll from 
the fractional differential equation describing a kicked 
system, is perhaps the best candidate to start a general 
investigation of the properties of maps which can be ob- 
tained from FDEs. 
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As it was shown in |15], maps that can be derived 
from FDEs are of the type of discrete maps with mem- 
ory. One-dimensional maps with memory, in which the 
present state of evolution depends on all past states, 
studied previously IH [I^ [ll, EUS Ee were not 
derived from differential equations. Most results were 
obtained for the generalizations of the logistic map. 

In the physical systems the transition from integer or- 
der time derivatives to fractional (of a lesser order) in- 
troduces additional damping and is similar in appear- 
ance to additional friction [a, ES]. Accordingly, in the 
case of the FSM we may expect transformation of the is- 
lands of stability and the accelerator mode islands into 
attractors (points, attracting trajectories, strange attrac- 
tors). Because the damping in systems with fractional 
derivatives is based on the internal causes different from 



the external forces of friction 11221 12311 . the correspond- 
ing attractors are also different from the attractors of 
the regular systems with friction and are called frac- 



tional attractors 

iBIBESIil 



Even in one-dimensional cases 



2111 most of the results were ob- 



tained numerically. An additional dimension makes the 
problem even more complex and most of the results in 
the present paper were obtained numerically. 

2. FSM, initial conditions 

The standard map in the form 

Pn+i ^ Pn-Ksinxn, x„+i =x„+p„+i (mod27r)(l) 
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can be derived from the differential equation 

oo 

'x + K s,\n(x) - »] = 0. 



(2) 



By replacing the second-order time derivative in eq. 
(|2]l with the Riemann-Liouville derivative off one ob- 
tains a fractional equation of the motion in the form 



CO 

D'^x+Ksin{x)^5i--n[ = 0, (1 < a < 2), (3) 

«=o 

oD^x(f) = D", oir^xit) = 



where 



1 



r(« - a) dt" Jo it - t) 



x{T)dT 



(« - 1 < ff < n), (4) 



D" - d"/dt", and qI^ is a fractional integral. The initial 
conditions for (O are 



ioDr'x)(0+)^pu 



(oD-J-^x)(0+) = ^7. 



(5) 



The Cauchy type problem ((Sj and (|5]l is equivalent to the 
Vol terra integral equation of the second kind 124125112611 



^ w sin 
" Jo ~ 



r(a) r(a - 1) 



Defining the momentum as 

Pit) = oz^r'^co, 



(6) 



(7) 



and performing integration in Q one can derive the 
equation for the FSM in the form (for the thorough 
derivation see 1I26I1 ') 



Pn+i = Pn - Ksinx,,, 
1 " 

Xn+1 = rpr— y,Pi+iV„in - i + l)+ 
via) ^ 



(8) 



1=0 



in + \ f-^, (mod 2n), 



r(a-i) 

where 



Vaim) = m""' -(m- If"'. 



(9) 



(10) 



Here it is assumed that T - \ and 1 < a <2. The form 
of eq. (|9]l which provides a more clear correspondence 



with the SM ia - 2) in the case b = is presented in 
Sec.EKeq. (EB)- 

The second initial condition in (jSj can be written as 



oZ)r'(0+) = lim olf-^xit) = 
r->0+ 

r'^^rWr (i<,<2),(ll) 
r-o+r(2-a)Jo (f-T)«-i 

which requires = in order to have a solution bounded 
at f = for a < 2. The assumption b - Q leads to 
the FSM equations which in the limiting case a - 2 
coincide with the equations for the standard map under 
the condition xq - 0. 

In this paper the FSM is taken in the form derived in 
ifTsIl which coincides with (O and ^ if b — 0. It is also 
assumed that xq = and the results can be compared to 
those obtained for the SM with xq = and arbitrary po. 
As a test, for the SM and for the FSM with a - 2 and 
the same initial conditions numerical calculations show 
that phase portraits look identical. 

System of equations ([8]l and (|9]l can be considered 
either in a cylindrical phase space (x mod 2n) or in un- 
bounded phase space. The second case is convenient to 
study transport. The trajectories in the second case are 
easily related to the first case. The FSM has no period- 
icity in p (the SM does) and cannot be considered on a 
torus. 



3. Stable fixed point 

The SM has stable fixed points at (0,27rn) for K < 
Kc - 4. It is easy to see that point (0, 0) is also a fixed 
point for the FSM. Direct computations using ([8]l and 
^ demonstrate that for the small initial values of po 
there is a clear transition from the convergence to the 
fixed point to divergence when the value of the param- 
eter K crosses the curve K - Kda) on Fig. la from 
smaller to larger values. 

The following system describes the evolution of tra- 
jectories near fixed point (0, 0) 

5p„+i = 6p„ - K6x„, (12) 

1 " 

Sxn+i = — — - } 6pi+iVain - i + 1). (13) 
r(a) ^ 



=0 



The solution can be found in the form 

n-l 



■>". = "» Z Ml;) (S)' <">«>• <"> 



The origin of the terms in parentheses, as well as the 
definition 
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will become clear in Sec.|5] Eqs. (fT2l i - ( fTSI l lead to the 
following iterative relationships 

n 

Xn+i,i = - '^in-m+lf'^Xm,i~u (0 < / < n),(17) 



(18) 



with the initial and boundary conditions 

= = (-1)", Pn+l,0^l, (19) 

Jc«+i,o = (n + l)""'- 

From (fTTT i and (fTST l it is clear that the series (fT4l i and ( fTsT i 
are alternating and it is natural to apply the Dirichlet's 
test to verify their convergence. This can be done by 
considering the totals 
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Figure 1 : Stability of the fixed point(0,0): a). The fixed point is stable 
below the curve K = Kc(a)', b). Values of 5c«j and /o<, obtained after 
20000 iterations of eqs. J221 and As o- — » 2 the values Soo and 
/oo increase rapidly. For a = 1.999, Soo ~ 276 and /„ » 552 after 
20000 iterations; c). An example of the typical evolution of Sco and 
loo over the first 200 iterations for 1 < o- < 2. This particular figure 
coiTesponds to = 1.8; d). Deviation of the values S„ and /„ from 
the values Soo « 2.04337 and /„ a 3.37416 for a = 1.8 during the 
first 20000 iterations (this type of behavior remains for 1 < a < 2); e). 
Evolution of trajectories with po = 1.5 + 0.0005;, < i < 200 for the 
case K = 3,a = 1.9. The line segments correspond to the nth iteration 
on the set of trajectories with close initial conditions. The evolution 
of the trajectories with smaller po is similar; f). 10^ iterations on both 
of two trajectories for K = 2, a = 1.4. The one at the bottom with 
Po = 0.3 is a fast converging trajectory. The upper trajectory with 
Po = 5.3 is an example of the ASCT in which piooooo ~ 0.042. 



They obey the following iterative rules 



n-l 



(21) 



S„^n"-' - — Vin-ir-'Si, = (22) 



(23) 



Computer simulations show that values of S„ and /„ 
converge to the values (-l)"^^Sc<j and (-l)"^'/oo de- 
picted on Fig. lb. Figs. Ic, Id show an example of the 
typical evolution of S„ and /„ over the first 20000 itera- 
tions. It means that the condition of convergence of 6p„ 
and Sx„ is 



ValK 

2r(a) 



< 1. 



(24) 



Numerical evaluation of the equality K - 2T{a)IVai 
perfectly reproduces the curve on Fig. la obtained by 
the direct computations of dS) and 

Because not only the stability problem (fTSl i and (fTsT l, 
but also the original map (jSj and (|9]l, contains convolu- 
tions, the use of generating functions 1,27] . which allows 
transformations of sums of products into products of 
sums, could be utilized in the investigation of the FSM 
and some other maps with memory. As an example, in 
the particular case of the stability problem (fT2l i and (fTsT l. 
the introduction of the generating functions 
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Pit) = Po- 



1 + Wait) 
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(l - Wait)) 



(25) 



(26) 



(27) 



(28) 



(29) 



Now the original problem is reduced to the problem of 
the asymptotic behavior at f = of the derivatives of 
the analytic functions Xit) and Pit), which is still quite 
complex and is not considered in this article. 

In the region of the parameter space where the fixed 
point is stable, the fixed point is surrounded by a finite 
basin of attraction, whose width W depends on the val- 
ues of K and a. For example, for K - 3 and a - 1.9 
the width of the basin of attraction is 1.6 < W < 1.7. 
Simulations of thousands of trajectories with pQ < 1.6 
performed by the authors, of which only 200 (with 
1.5 < pd < 1.6) are presented in Fig. le, show only 
converging trajectories, whereas among 200 trajecto- 
ries with 1.6 < pq < 1.7 in Fig 2a there are trajecto- 
ries converging to the fixed point as well as some tra- 
jectories converging to attracting slow diverging trajec- 
tories (ASDT), whose properties will be discussed in 
the following section. Trajectories in Fig. le converge 
very rapidly. In the case K - 2 and a - 1.4 in addi- 
tion to the trajectories which converge rapidly and AS- 
DTs there exist attracting slow converging trajectories 
(ASCT) (Fig. If). 
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Figure 2: Phase space with ASDTs: a). The same values of parame- 
ters as in Fig le but po = 1.6 + 0.0005!; b). 200 iterations on trajec- 
tories with po = 4 + 0.02i, < ! < 500 for the case K = 2, a = 1.9. 
Trajectories converging to the fixed point, ASDTs with x = 0, and 
period 4 attracting trajectories are present; c). 2000 iterations on tra- 
jectories with Po = 2 + 0.04;, < r < 50 for the case K = 0.6, 
a = 1.9. Trajectories converging to the fixed point, period 2 and 3 
attracting trajectories are present; d).The same values of parameters 
as in Fig le but po = 5 + 0.005!. 



4. Attracting slow diverging trajectories (ASDT) 

As it can be seen from Fig 2a, the phase portrait on a 
cylinder of the FSM with K - and a - 1.9 contains 
only one fixed point and ASDTs approximately equally 
spaced along the p-axis. This result corresponds to the 
fact that the standard map with K - 2> has only one 
central island. More complex structure of the standard 
map's phase space for smaller values of K (for exam- 
ple for K - 2 and K - 0.6) can explain more complex 
structure of the FSM's phase space, where periodic at- 
tracting trajectories with periods T - A (Fig. 2b), T -2, 
and r = 3 (Fig. 2c) are present. 

Each ASDT has its own basin of attraction (see 
Fig. 2d). Between those basins two initially close tra- 
jectories at first diverge, but then converge to the same 
or different fixed point or ASDT. 

Numerical evaluation shows that for ASDTs which 
converge to trajectories along the p-axis (jc — > jc/,,,, - 
0) in the area of stability (which is the same as for the 
stability of the fixed point) the following holds (for large 
n see Fig. 3a) 



Pn 



Cn 



2-a 



(30) 



The constant C can be easily evaluated for 1.8 < a < 2. 
Consider an ASDT with jc/,,,, - Q, T - \, and 2nM, 
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Figure 3: Evaluation of the behavior of the ASDTs: a). Momenta for 
two ASDTs with x„ k 2nn in the unbounded space (in this example 
K = 2). The solid line is related to a trajectory with a = 1.9 and its 
slope is 0. 1 . The dashed line corresponds to a trajectory with a = \.5 
and its slope is 0.5; b). Deviation of momenta from the asymptotic 
formula for two ASDTs with x,, x 2nn in the unbounded space, a = 
1.9, and K = 2. The dashed line has po = 1 and the solid one po = 6; 
c). Relative deviation of the momenta for the trajectories in Fig 3b 
from the asymptotic formula; d). Deviation of the x-coordinates for 
the trajectories in Fig 3b from the asymptotic formula. 



where M is an integer, constant step in x in the un- 
bounded space. Then Eq. (|9]l with b - Q gives 



1 " 

• x„ = — — y^ {pk+\ - PkWain - ^ + 1) (31) 

1 (IT) ^—^ 

■Vain + 1), 



k=l 

Pi 

'T(a) 



For large n the last term is small (~ n" ^) and the fol- 
lowing holds 

n 

J](PM - Pk)ycin - + 1) = 2nMY{a). (32) 

k=\ 

With the assumption p„ ~ n^'" it can be shown that 
for values of a > 1.8 considered the terms in the last 
sum with large k are small and in the series representa- 
tion of Va(n - k + 1) it is possible to keep only terms of 
the highest order in k/n. Thus, (l32T i leads to the approx- 
imations 



p„ ~ po + 



a- 1 
InMil - a)r(a) 



(33) 



(34) 



In the case K - 2, a - 1.9 Figs. 3b-3d show two trajec- 
tories with M = 1 (initial momenta po = 6 and po = 7) 



approaching an ASDT: the deviation from the asymp- 
totic ( l33T l and ( l34l i and the relative difference with re- 
spect to ( l33T l. 



5. Period 2 stable trajectory 

The SM has two stable points of the period T - 2 
trajectory for A < K <2n with the property 



Pn+\ 



Pm X„+i 



(35) 



The same points persist in the numerical experiments 
for the FSM (Fig 4a). These points are attracting most 
of the trajectories with small po. Assuming the exis- 
tence of a r = 2 attracting trajectory, it is possible to 
calculate the coordinates of its attracting points {xi,pi) 
and {-xi, -pi). In this case from ([8]) and ^ 



P'^ 2 



K 



2r(a) 



lixdJ^i-lf^'Vaik) 



(36) 



(37) 



k=l 



Finally, the equation for xi takes the form 
K 



Xl 



where 



2r(a) 



Vai sin(x/). 



(38) 



(39) 



k=l 



and can be easily calculated numerically. From ( |38] ) the 
condition of the existence of T = 2 trajectory 



2r(a) 

Vai ' 



(40) 



is exactly opposite to ( l24l i. It is satisfied above the curve 
K = Kc(a) on the Fig. la. For a = 2 (l40b produces the 
well-known condition K > 4 for the SM. The results 
of calculations of the xi and pi for the cases K - 4.5, 
1 < Q- < 2 presented in Fig. 4b-d perfectly coincide 
with the results of the direct computations of ([8]l and (|9]l 
with b = 0. After 1000 iterations presented in Figs. 4e,f 
the values of deviations \pn - pi\ and \xn - xi\ are less 
than 10"^ 
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Figure 4: Period 2 stable trajectory: a). An example of r = 2 attractor 
for K = 4.5, a = 1.9. One trajectory with xq = 0, po = 0.513; b). pi 
of xi for the case of *r = 4.5; c). pi of a for the case of ^ = 4.5; d). xi 
of a for the case of ^ = 4.5; e). p„ - pi for the trajectory in Fig. 4a. 
After 1000 iterations \p„ - pi\ < 10 '; f). x„ - xi for the trajectory in 
Fig. 4a. After 1000 iterations \x„ - xi\ < lO"''. 
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Figure 5: Cascade of bifurcation type trajectories: a). 120000 itera- 
tions on a single trajectory with K = 4.5, a = 1.65, po = 0.3. The 
trajectory occasionally sticks to a CBTT but then always recovers into 
the chaotic sea; b). 100000 iterations on a trajectory with K = 3.5, 
a = 1 . 1 , po = 20. The trajectory very fast turns into a CBTT which 
slowly converges to a fractal type area. 
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Figure 6: Examples of phase space for K > 1w. a). An attracting 
ballistic trajectory with K = 6.908745, a = 1.999, po = 0.7; b). A 
chaotic trajectory for K = 6.908745, a = 1.9. 

6. Cascade of bifurcations type trajectories (CBTT) 

Period 2 stable trajectories have limited basins of 
attraction. Trajectories that don't fall into those ar- 
eas reveal a diverse variety of properties, from period 
two slow attracting trajectories to fractal type attractors 
and cascade of bifurcations type trajectories (CBTT). 
Fig. 5a presents a single chaotic trajectory which sticks 
to the areas similar to the cascade of bifurcations which 
are well-known for the logistic map. In Fig. 5b a sin- 
gle trajectory falls very rapidly into one of the attracting 
CBTTs. Because the bifurcation diagram of the logistic 
map has fractal properties (see for example Chapter 2 in 
1280), it is expected that the structure to which this tra- 
jectory slowly converges also possesses fractal features. 

The properties of this type of attractors, as well as the 
properties different types of observed during computer 
simulations chaotic, attracting, and ballistic trajectories 
for K > 2n (see Fig 6) will be considered in the subse- 
quent article. 



by a relatively small number of computations and this 
is understandable, taking into account all the difficulties 
of performing numerical simulations for the equations 
with fractional derivatives. 

The use of the FSM, which is equivalent to the origi- 
nal FDE, allows performing thousands of runs of simu- 
lations of the kicked fractional system with two param- 
eters: K and a. The FSM also allows making some an- 
alytic deductions and revealing some properties of the 
fractional attractors which were not reported before: 

a) . The stability of the fixed point (0,0) of the FSM 
is different not only from the stability of the fixed point 
in the domain of the regular motion (zero Lyapunov ex- 
ponent) of the SM, but also from the stability of fixed 
attracting points of the regular (not fractional) dissipa- 
tive systems like, for example, the dissipative standard 
map (Zaslavsky map) foil . The difference is in the way 
in which trajectories approach the attracting point. In 
the FSM this way depends on the initial conditions. For 
example, in Fig. If there are two trajectories approach- 
ing the same fixed point: one is fast spiraling into the 
attractor and the other is slowly converging. 

b) . Stable period 2 attracting trajectories exist only 
in the asymptotic sense - they do not represent any real 
periodic solutions. If the initial condition is chosen in 
a period two stable attracting point, this trajectory will 
immediately jump out of this point and where it will end 
depends on the values of K and a. 

c) . All the FSM attractors exist in the sense that 
there are trajectories which converge into those attrac- 
tors. But if an initial condition is taken on any of the 
attracting trajectories (except for the fixed point), they 
will most likely not evolve along the same trajectory. 



7. Fractional attractors and their stability 

The problems of existence and stability of the frac- 
tional attractors for the systems described by the FDEs 
were addressed in a few recent papers. It was noticed in 
II22II that the properties of the fractional chaotic attrac- 
tors are different from the properties of the "regular" 
chaotic attractors and may have some pseudochaotic 
features. The problem of existence of multi-scroll frac- 
tional chaotic attractors was considered in fl^. The 
problem of stability of the stationary solutions (fixed 
points for ODEs) of systems described by the fractional 
ODEs and PDEs was considered in f30l |3]J, |32|]. In 
the above mentioned articles the equations contained the 
Caputo fractional derivatives, whereas in the present ar- 
ticle the Riemann-Liouville fractional derivative is used. 
This fact does not allow a direct comparison of the re- 
sults. The results Q IH [l^, IH [ll were supported 



8. Conclusion 

In this article properties of the phase space of the 
FSM were investigated. It was shown that islands of 
regular motion of the SM in the FSM turn into at- 
tractors (points, attracting trajectories, and fractal-like 
structures). Properties of the attracting fixed points, pe- 
riod two trajectories, ASCTs, and ASDTs were consid- 
ered. This consideration allows the description of the 
evolution of the dynamical variable x of the original 
fractional dynamical system, a system described by the 
FDE reducible to the FSM. Physical interpretation of 
the momentum, defined through a fractional derivative 
from the variable x, is unclear 

The explanation of the CBTTs, which are interest- 
ing phenomena, requires further detailed investigation. 
Chaotic trajectories that spend some time near CBTTs, 



which can be called sticky attractors in analogy to sticky 
islands of the SM, are good candidates for the investi- 
gation of anomalous diffusion. Transport was not con- 
sidered in this article. How general the properties of 
the phase space of the FSM are will become clear after 
further investigations of different fractional maps, maps 
with memory which can be derived from the FDEs, and 
particular those suggested in [lE], will be conducted. 
The fact that so many physical systems can be reduced 
to studying of the SM gives a hope that those physical 
systems which can be reduced to studying the FSM will 
be found. 
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